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1. Introduction 


A family of open sets with smooth boundary {Et}o<t<T in is said to 
move according to volume-preserving mean-curvature flow if the motion law, 
expressed as an evolution equation for the boundaries dEt, takes the form 


V = {H)-H on dEt, (1.1) 

for all t G [0,T]. Here, at any point x on dEt, v{x) denotes the velocity 
component normal to the boundary, in the direction of the outer normal, H(x) 
is the scalar mean curvature (with the sign convention that H is positive for 
balls, see the next section), and the brackets (•) denote the average of a quantity 
over the boundary of Et- 


It is immediately verified that the volume of the sets Et (i-e., its n-dimensional 
Lebesgue measure, denoted by \Et\) is indeed preserved under the smooth flow 
(1.1) because 


dt 


\Et\ = 


I 

JdEt 


V dP 


n-l (1-1) 


= ' 0 . 


And thus, upon rescaling variables, we may assume that \Et\ = |Eo| = 1 for 
any t G [0,T]. Moreover, the perimeter of the sets Et is decreasing because 


^7r-^{dEt) = f vHdn^-^ ^ = ^ 

dt ^ Jqe, 



< 0 . 


During a typical evolution, a volume-preserving mean-curvature flow exhibits 
singularities of different kinds, even in the case of smooth initial data. These 
singularities correspond to changes in the topology of the configuration and 
include shrinkage of islands to points and disappearance, collisions and merging 
of neighboring islands, pinch-offs etc... If the topology changes, the boundary 
of the evolving set looses regularity and, as a consequence, the formulation 
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(1.1) of the evolution law is inadequate. The goal of the present work is the 
construction of a notion of a weak solution to the volume-preserving mean- 
curvature flow that is global in time and thus overcomes these singular moments. 

Several solutions to volume-preserving mean-curvature flow have been pro¬ 
posed in the literature: existence and uniqueness of a global in time smooth 
solution and its convergence to a sphere is shown in [17, 22] for smooth convex 
initial data and in [15, 24] for initial data close to a sphere (for further re¬ 
lated results see [4, 7, 8] and the references therein). In principle, these results 
also yield local in time existence and uniqueness of smooth solutions. In [30] 
and [31], the authors consider level-set and diffusion-generated solutions for the 
purpose of numerical studies. In [5], Bellettini, Caselles, Chambolle and No- 
vaga construct solutions for the volume-preserving anisotropic mean-curvature 
(or anisotropic variant of (1.1)) for convex sets using a method similar to ours 
(that will be outlined below). In this paper, it is also shown that the solution 
(the so-called flat flow) is unique and coincides with regular solutions when the 
latter are defined. We also mention a mean field approximation approach to 
the volume preserving mean-curvature flow as developed in [6, 12, 1]. 

It is well-known that volume-preserving mean-curvature can be (formally) 
interpreted as the L^-gradient flow of the perimeter functional for configura¬ 
tions with a fixed volume, see, e.g., [29, Sec. 2]. This gradient flow structure, 
however, is for the purpose of well-posedness results impracticable, since the 
(geodesic) distance is degenerate in the sense that two well separated conhgura- 
tions may have zero distance [28]. In the present manuscript, we follow the 
method proposed independently by Almgren, Taylor & Wang [2] and Luckhaus 
&: Sturzenhecker [26] in the study of (forced) mean-curvature flows to bypass 
this difficulty. The authors consider an implicit time-discretization of the flow, 
which comes as a gradient flow of the perimeter functional with respect to a 
new non-degenerate distance function that approximates the distance. The 
limiting time-continuous flow constructed with this method is usually referred 
to as the flat flow. The main difference between the present work and [2, 26] 
relies on the non-locality of the volume-preserving mean-curvature flow. As 
an immediate consequence, there is no maximum-principle available for (1.1). 
Also related to this aspect, there is the problem of proving the consistency of 
the scheme, i.e. the coincidence of the flat volume-preserving mean-curvature 
flow with the smooth one when the latter exists. Under some assumption on 
the Lagrange multiplier of the flow, the consistency can be inferred from the 
arguments in [2] , but we do not know if these conditions are generally satished 
and we do not discuss further the problem of the consistency in this paper. A 
more detailed discussion on the different features of the flows in [2, 26] and the 
one considered in the present manuscript will follow in Section 3 below. 

We conclude this subsection with a short discussion on the background of this 
evolution. Volume-preserving mean-curvature flow can be considered as a sim¬ 
plified model for attachment-limited kinetics, and as such it plays an important 
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role in the study of solidification processes, where solid islands grow in an under¬ 
cooled liquid of the same substance. In such situations, solid particles melt at 
high-curvature regions and simultaneously precipitate at low-curvature regions, 
while the total mass of the solid remains essentially constant [33, 11, 32]. In this 
way, the total surface area of solid islands is decreasing, and thus, this process 
leads to the growth of larger islands at the expense of smaller ones: a phe¬ 
nomenon called coarsening [29]. In general, solidification processes are mathe¬ 
matically often modelled by Mullins-Sekerka equations (or a Stefan problem), 
where the Gibbs-Thompson relation is modified by a kinetic drift term [25, 21], 
and their phase-field counterparts respectively [9]. This model allows for both 
attachment kinetics (kinetic drift) and bulk diffusion (Mullins-Sekerka/Stefan). 
It turns out that attachment kinetics is the relevant mass transport mechanism 
in earlier stages of the evolution while bulk diffusion predominates the later 
stages [13]. In a certain sense, volume-preserving mean-curvature flow natu¬ 
rally arises as the singular limit of this more general solidification model in the 
regime of vanishing bulk diffusion. More recently, variants of volume-preserving 
mean curvature flow were also applied in the context of shape recovery in image 
processing [10]. 

The article is organized as follows: in Section 2 we fix the notation and state 
the main results of the paper, which are then proved in Sections 3 and 4 and are 
the existence of flat volume-preserving mean-curvature flows and the existence 
of distributional solutions, respectively. 
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2. Statements of the main results 

2.1. Notation. For any Lebesgue measurable set E C K”, we denote by \E\ 
the n-dimensional Lebesgue measure of E and by xe the characteristic function 
of E, i.e. xe{x) = 0 if X ^ E and xe{x) = 1 li x ^ E. The perimeter of E in 
an open set C ffi"' is defined as 

Per (FI, 11) := sup < / divipdx: y? G (7^(11; M") with sup |:^| < 1 i , 

[Je n } 

and we write Per(Fl) := Per(i?, M”). If the latter quantity is hnite, we will 
call E a set of finite perimeter. In the case that E is an open set with dE of 
class C^, we simply have Per(Fl, Q) = T-E~^{dE n VL) and Per(Fl) = T-E~^{dE). 
The reduced boundary of a set of finite perimeter E is denoted by d*E, cp. 
[16, Sec. 5.7], and for the unit outer normal to E we write i'e- The tangential 
divergence of a vector field 'L G C'^(M”';M”) with respect to dE is defined by 
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divg^;'!' := div^ — ue ■ We say that a set of finite perimeter E has a 

(generalized) mean-curvature He S {d*E, provided that 


/ ' 

Ja*E 


diva^;'!' (HE ^ = 


'h • ueHe <HL 


n—l 


for all 'h G C, 


1 (Tan. Tn)n\ 


ld*E 


( 2 . 1 ) 


Observe that with this sign convention it is //br = ^ 


R ■ 


We write sdi;’ for the signed distance from a Lebesgue measurable set E, with 
the convention that sdj? is negative inside E, i.e., 


sdFix) 


dist(x,F) for X € 

— dist(x, F‘^) for x G F, 


Here, F‘^ := M” \ F denotes the complement set of F, and the distance from a 
set F is by definition 

dist(x,F) = inf \x — y\. 

y&F 

We will sometimes use the notation di? = Isd^l. 

By [t] we denote the integer part of a real number t, that is the biggest integer 
m such that m < t. 


Finally, we denote by Cn,Cn positive constants that depend on the space 
dimension only. Moreover, Cn,o, Cn,o, and Cn,o,T are constants that may addi¬ 
tionally depend on the initial data or the time T > 0. During the computations, 
the value of these constants may change from line to line. However, for the sake 
of clarity we need to keep track of the dimensional constant in Proposition 3.2.1, 
therefore we make an exception to the above convention and denote it by 7 ^. 
The volume of the n-dimensional unit ball will be denoted by ujn, and thus its 
surface area is mon- 


2.2. Approximate solutions. In this paper we introduce a notion of global 
flat solution to the volume-preserving mean-curvature flow which is based on 
the implicit time discretization of (1.1) in the spirit of Almgren, Taylor & Wang 
[2] and Luckhaus & Sturzenhecker [26]. That is, we consider a time-discrete 
gradient flow for the perimeter functional. For this purpose, we define 

Eh{E,F) :=Per(F) + i^sdFdx + ^||F| -1|, 

for any two sets of finite perimeter E and F in Here, h is a positive small 
number that plays the role of the time step of approximate solutions. The 
second term in the above functional approximates the degenerate geodesic 
distance on the configuration space of hypersurfaces. The functional differs 
from the one considered in the original papers [2, 26] only in the last term: a 
weak penalization that favors unit-volume of minimizing sets. 
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2 . 2 . 1 . Definition. Let Eq be a set of finite perimeter with \Eq\ = 1, and 

h> 0. Let be a sequence of sets defined iteratively by 

E^q^ = Eq and E^^j^ G argmin | for k > 1. 

_BCR" ) J 

We furthermore define 

E^'^ := E^^jI for any t G [kh, {k + l)h), 

and call }t>o an approximate flat solution to the volume-preserving mean- 
curvature flow with initial datum Eq . 

The existence of minimizers E^ and thus the existence of an approximate 
flat solution is guaranteed by Lemma 3.1.1 below. Incorporating the volume 
constraint in a soft way into the energy functional rather than imposing a 
hard constraint on the admissible sets has the advantage that we are free to 
chose arbitrary competitors, most notably in the derivation of density estimates. 
Thanks to the penalizing factor 1/^/h, the constraint becomes active in the limit 
h fO. Even more can be shown; the number of time steps in which approximate 
solutions violate the volume constraint \E^^'^\ = 1 can be bounded uniformly 
in h, cf. Corollary 3.4.5 below. A similar functional including a soft volume 
constraint was recently considered by Goldman &: Novaga [20] in the study of 
a prescribed curvature problem. 

2.3. Main results. We can now state our main results. The first one is a 
convergence result for approximate solutions. 

2 . 3 . 1 . Theorem (Existence of flat flows). Let Eq be a bounded set of finite 

perimeter with jEoj = 1 and, for any /i > 0, let be an approximate flat 

solution to the volume-preserving mean-curvature flow with initial datum Eq. 
Then, there exists a family of sets of finite perimeter {Et}t>Q and a subsequence 
/ifc i 0 such that 

—>• 0 for a.e. t G [0,-|-oo), 

and, for every 0 < s < 

\Et\ = 1 , 

\EtAEs\ < Cn,o\t - s\^/‘^, 

Per(Et) < Per(Es). 

Our next statement is the existence of a distributional solution in the sense of 
Luckhaus &: Sturzenhecker [26] to the volume-preserving mean-curvature flow 
under the hypothesis that the perimeters of the approximate solutions converge 
to the perimeter of the limiting solutions identified in the previous theorem. 
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2.3.2. Theorem (Existence of distributional solutions). Suppose that n <7. 
Let (/ifc)fceN {Et}t>o be as in Theorem 2.3.1. For any T > 0, if 


lim 

k^oo 


nT 

/ Per(Ef''^) 

Jo 


dt = 


Per^Et) dt, 


( 2 . 2 ) 


then {Et}o<t<T is a distributional solution to the volume-preserving mean- 
curvature flow with initial datum Eq in the following sense: 

(1) for almost every t € [0,T) the set Et has (generalized) mean curvature 
in the sense of (2.1) satisfying 

rT, 


[ [ \He-. 

Jo Jd*Et 


< +oo; 


(2.3) 


(2) there exists v : M"’ x (0,T) — )• R with v{-,t)\g*Et G L(^{d*Et,dFE ^) for 
a.e. t € (0,T) and Jq fg*Et v^dFL^~^dt < +oo such that 


-[ [ v(fd'H^-^dt= f [ iHEtf-Xf))d'H'"-^dt, 

Jo Jd*Et Jo Jd*Et 

n dtfdxdt + f 4>{0, ■) dx = — f f 
J En J 0 J d 


dt4> dxdt 

lo JEt 

for every f £ (^^([OjOo) x R”), where 
1 


vfdFE ^dt, 


0 Jd*Et 


m ■■= 


n^-\d*Et) Ja^E, 


HE,dH 


n—1 


for a.e. t £ [0,T). 


(2.4) 

(2.5) 

( 2 . 6 ) 


In the second part of the theorem, Lq is the set of all functions with zero 
mean. 

Note that (2.4) is a weak formulation of (1.1), while (2.5) establishes the link 
between v and the velocity of the boundaries of Et. It is straightforward to 
check that smooth solutions of (1.1) satisfy (2.4) and (2.5). 

The restriction on the dimension n < 7 is technical and is needed in the 
proof of Corollary 4.2.1 where the Bernstein theorem for minimal surfaces is 
exploited. 


3. Flat volume-preserving mean-curvature flows 

In this section we prove the first main result in Theorem 2.3.1. We follow 
quite closely Luckhaus &: Sturzenhecker [26], providing all the details for the 
readers’ convenience. 


3.1. Existence of approximate solutions. We start remarking that 


sd_F dx = 


/ d-F dx — 
Jeaf Jf 


dF dx. 


(3.1) 


The existence of the approximate flat solutions is guaranteed by the following 
lemma. 
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3 . 1 . 1 . Lemma. Let F C M” be a bounded set of finite perimeter. For every 
h > 0, there exists a minimizer E of , F) and, moreover, E satisfies the 
discrete dissipation inequality 

Per(E) + 11 dFdx + ^\\E\-l\ < Per(F) + ^ ||F| - 1| . (3.2) 

hjEAF \h Vh 

Proof. Since E is an admissible competitor, we obtain by (3.1) that 
0 < mfFh{E,F) + — I dp dx < Per(F)-|—^||F| —1| < oo. (3.3) 

E h Jf y/h 

Let denote a minimizing sequence oi Fh{ - ,F). Without loss of gen¬ 

erality we may assume that Ey dC. Bp for a suitable i? > 0. Since {x_E,,}jyeN is 
bounded in BV{Bf), there exists a subsequence (not relabeled) that converges 
weakly to a function x iii BV{Be), and thus strongly in L^{Bf). In partic¬ 
ular, X is the characteristic function of some set of finite perimeter E. Since 
Xe sdi? dx is continuous and the perimeter is lower semi-continuous with 
respect to convergence, it follows that 

Fh{E,F) < liminf Fh{Ey, E) = inf Fh{E,F). 
i/foo EcK" 

Therefore, E minimizes Fh{- ,F) and (3.2) follows from (3.3). □ 

By standard results on minimal surfaces (see [27]), it holds that the mini- 
mizers E of Fh{ ■ , F) can be chosen to be closed subsets with dE of class 
up to a (relatively closed) singular set of dimension at most n — 7. Using the 
Euler-Lagrange equation for Fh{ ■ ,F), one can also show that the regular part 
of the boundary dE is actually (cp. Lemma 3.4.2). 

3.2. and L^-estimates. Our next statement gives a uniform bound on 
the distance between the boundary of the minimizing set and the boundary of 
the reference set. 

3 . 2 . 1 . Proposition (L°°-estimate). There exists a dimensional constant 
7„ > 0 with the following property. Let F C M"" be a bounded set of finite 
perimeter and let E be a minimizer of Fh{ ■ ,F). Then, 

sup dE < JnVh. (3-4) 

EAF 

Proof. The proof of this proposition is based on the density estimates for 
one-side minimizing set which for readers’ convenience we prove in the Appen¬ 
dix A. We claim indeed that the statement holds with 

r dnujn\ 

7 „ = max I 3,-I, 

I Cn J 

where Cn is the dimensional constant in Lemma A.0.1. The argument is by 
contradiction. Let c > max|3, | and let xq G EAE contradict (3.4) with 
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7 „ replaced by c. Without loss of generality, we can assume that xq & F \ E: 
the other case is at all analogous. We then have that 

sdi?(xo) < —cVh. (3-5) 

Then any ball Br{xo) of radius r < is contained in F. By the minimality 
of E, we have Fh{E, E) < Fh{E D Br{xQ), E), and thus 


If 1 

Per(£i) < Per(ill U i?r(a;o)) +-r / sdp dx-\ —^|i?r.(xo) \ (3-6) 

d JBr{xo)\E Vh 


We use (3.5) and r < to infer that 


1 

h 


/ sdp dx < - -j= 

'Br{xo)\E 2Vh 


Br{xo) \ E\. 


(3.7) 


Then (3.6) and (3.7) yield 

Per(.E) < PeiiE U Br{xo)) - - l) \Br{xo) \ E\. (3.8) 

By assumption c > 3 and we can apply Lemma A.0.1 with = 0 and obtain 

cVh 


\Brixo) \ Lll > CnV^ for a.e. 0 < r < 


(3.9) 


On the other hand, from (3.8) we deduce also that for a.e. 0 < r < 

< n^-^{dBr{xo) \E)< (3.10) 


Combining (3.9) and (3.10), we get that 

Cnr"'<\Br{xo)\E\ <nu}n Vhr'^~^, 

for almost all 0 < r < which gives the desired contradiction to the choice 
of c as soon as r t • □ 

The following density estimates are now an immediate consequence. 


3.2.2. Corollary. Let E C M"" be a bounded set of finite perimeter and let 
E be a minimizer of Fh{-, E). Then, for every r £ (0,7„\//i) and for every 
xq £ dE, it holds 

mm{\Br{xo)\E\,\E n Brixo)\} ^Cnr"^, ( 3 . 11 ) 

cnr^-^ < F&{E,Br{xo)) < Cnr^-\ ( 3 . 12 ) 

Proof. Since 77 is a minimizer of Fh{ ■ ,F), for any xq £ dE, it holds that 
Fh{E, E) < Fh{E \ Br{xo), E), which implies 

If 1 

Per(£') + — / sdi? dx < Per(i? \ ilj.(a;o)) H —-^\E Ci Br{xo)\. 

djEnBrixo) Vh 
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Estimating the second term via Proposition 3.2.1, we obtain 

n 

Per(E) < Per(£' \ Br{xo)) H —^\E n Br{xo)\. (3.13) 

v/i 

A similar analysis shows that 

Per(E) < Per(E U Br{xo)) + %|5,(xo) \ E\. 

Vh 

Therefore, by Lemma A.0.1 we deduce (by possibly redefining Cn) 

min I IE n Br{xo)\, \ Br{xo) \ E|} > c„ r” V 0 < r < jn'/h- 

The first inequality in (3.12) is now an immediate consequence of the relative 
isoperimetric inequality (cf. [19, Cor. 1.29]). For the second inequality, we 
rewrite (3.13) as 

n^-^dE n Brixo)) < n^-^dBrixo) n E) + %|E n E^(xo)|. 

y/h 

Since r < the upper bound is obvious. □ 

Next we prove an estimate on the volume of the symmetric difference of two 
consecutive sets of the approximate solutions. 

3.2.3. Proposition (L^-estimate). Let F C M"' be a bounded set of finite 
perimeter and let E he a minimizer of Fh{ ■ ,E). Then, 

\EAF\ < Cn(iPer{E) + l [ dpdx] (3.14) 

V « Jeaf J 

Proof. In order to estimate EAE, we split it into two parts: 

|EAE1 < |{x e EAE : dpix) < l}\ + \{x G EAE : dF(x) > £}\. 

The second term is easily estimated by 

|{x G EAE : di7’(x) > £}\ < - f dF{x)dx. 

^ Jeaf 

To estimate the first term, we use a simple covering argument to find a col¬ 
lection of disjoint balls {Bi{xi)}i£j with Xi £ d*E and I C N a finite set such 
that d*E C {JiQiB 2 e{xi). Note that by (3.11) and the relative isoperimetric 
inequality (cf. [19, Cor. 1.29]) we have for every i £ I 

iBsfixi)] < C„min{|E n E£(xj)l, |E£(a;j) \ E|} 
<Cnin^-\d*EnBe{xi)). 
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Note finally that the set {x G EAF : di7’(x) < is covered by 
Summing over i and the choice of the balls {Bi{xi)}i^i yields 

|{x G EAF : dpix) < £}| < ^ \B3i{xi)\ 

i&I 

iei 

<CniPev{E). □ 


3.3. Holder continuity in time. As an immediate consequence of the dis¬ 
crete dissipation inequality (3.2), we remark that 


Per(£', 


(h). 








dx H-= I E\ 




t — h 


y/h' 


<Per(E™) + ^||E 


(h) 


1| 

l| V t G [/i,-|-oo), (3.15) 


and, by iterating (3.15) and using iP^ol = 1, 


Per(A;f^) < Per(.Bo) V t > 0, 

J' I I T7i(^) 1 


y/h 

fl 


— 1 < Per(£'o) V t > 0, 


^t-h 




dx < Per(£'o), 


(3.16) 

(3.17) 

(3.18) 


for every T > h. Similarly, using Proposition 3.2.3 with i = h < jnVh, we also 
get 



W A TpW 


\Ey^’AE_ 


t—m 



<Cn{T + l)Pev{Eo). 


(3.19) 


3.3.1. Proposition (C^/^ regularity in time). Let h < 1 and let {Et}t>Q be 
an approximate flat flow. Then it holds 


_ ^|l/2 V 0 < t < S < +00. 

Proof. Clearly it is enough to consider the case s — t > h. Let j G N and 
k G N\{0} be such that t G [jh, {j + l)h) and s G [{j + k)h, {j + k + l)h). Then, 
we can use Proposition 3.2.3 with i = (note that i < jn'/h by the 
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assumption s — t > h) and (3.16), and estimate in the following way: 

k 


Er’AEfi < ^ ie;;VjAe, 




m=l 


{j+m)h 


|t-s|l/2 
m=l ' ' 






m=l 


L 




(j+m)h {j+m—l)h 


d (/i) dx. 


By using (3.15) we estimate the sum above by 


<CnYl 


h 


|i-s|V2 


Per(p;o) 


+ 


c„ E l‘ - »I‘''“(p«'(-e;‘1„-„J - 


m=l 

k 


+ Vh 


<Cn 


m=l 

kh 




{j+m-l)h 


- 1 - \\E 


{h) 

{j+m)h 


- 1 


|t-s|V2 


Per(.Eo) + Cn\t - s\^/\Vei{Et) - Per(P;,)) 


+ Cr. 


|t-s|V2 

y/h 


I Ply-1 — 1 — I i?., I — 1 


(3.20) 

Therefore, by (3.16) and (3.17), we get 

<Cn\t- s|'/2per(^o), (3.21) 

where we used kh < \t — s\ + h < 2|t — s|, thus concluding the proof of the 
proposition. □ 

3.4. First variations and first consequences. We now introduce the time- 
discrete normal velocity: for all t > 0 and x € M”, we set 


0 


^sd (/i) (x) for t G [h,+oo), 

^t—h 

for t G [0, h). 


3.4.1. Lemma (L^.bound on the velocity). Let he an approximate 

flat flow. Then it holds 


roo r 

Jo JdEr> 


(3.22) 


Proof. We fix t G [h, -|-oo) and consider for every £ G Z with 2^ < ’y^^/Vh the 
sets 

K{£) := {x G : 2^ < |x('*)(t,x)| < 2^+^}, 
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so that M” = It follows from 2^ < 7„\//i/2 and from Corollary 3.2.2 

that, for every x € dE[^\ 


nn-l^dE['^) < Cn(2^-^h) 


n—1 


(3.23) 

(3.24) 


Using 2^ ^ < \v^^'^{t,y)\ <4-2^ ^ for all y G B 2 t-ih{x) with x G dE[^^ n K{€}, 
we obtain for every x G dEf^^ n K{t} 


I 

[ C „( 2^-^)2 ( 2 ^- 1 / i ) 


\v^^^\dy Cn2^-^ (2^-^hy , 


Hence, combining these two estimates, we have 


[ (v^^^fdn^-^ < % / 

« iSy.l 




v^^'^ldy. 


Now, by a simple application of Besicovitch’s covering theorem [16, Ch. 1.5.2] 
to {i? 2 ^-i/j(a:) : x G dE^ n K{£)}, we obtain 

Cr, 


\v^^'^\dx. (3.25) 


[ U'^ydn^-^< ^ , 

Finally, summing up over £ G Z with 2^ < in (3.25) yields 

f (n(^))2d^"-i < ^2 f Iv^'^^ldx. 

We now show how the above estimate implies (3.22). In view of (3.15) we 
have 


JdEi’^^ 


2 dn^-^ 


< 




h \ 


- 1| - Per(£;f) - -^\\Ey\ - 1 


Integrating in time and using iFioj = 1, we obtain 

rT 


[ f dU^-^dt 

Jo JdE^^'^ 

< Cn(pej:{Ey)-Pev{Ey) 


E^J!^h - 1 


< C'nPer(£:o), 

from which, by a simple limit for T —)• + 00 , (3.22) follows. 


□ 
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We now derive the Euler-Lagrange equations which constitute the weak mo¬ 
tion law for the time-discrete evolution. 


3.4.2. Lemma (Euler-Lagrange equations). For every t G [/i,-|-oo) and 4' G 

1 


it holds 


[ fdiv -I-(i.) • 4') d'H"' ^ f v (h) 

JqeW \ dEl El > J El 


'i’d'H 


n—1 


IdE., 
where 


:= 




(3.26) 


Moreover, if \e[^\ / 1, then it also holds X^^\t) = isgn(l — \eI^^\). 


Wi 


Vh' 


As we shall see in the proof below, the constants X^^\t) defined in (3.27) 
are Lagrange multipliers corresponding to the volume constraint, whenever it is 
active. Since this constraint is satisfied up to a finite number of times (uniformly 
in h) by Corollary 3.4.5 below, by a slight abuse of terminology, we call X^^\t) 
a Lagrange multiplier, even if the volume constraint is not active. 


Proof. If \E^\ ^ 1, it is very simple to compute the variations of 
Fh{-,E^^) along the vector field 'k G C^(M”;M”) and see that they are given 
by (3.26) with X^^\t) = ;^sgn(l — \Ej:^^\). In the case \E^^'^\ = 1, we have 


G argmin < Per(E) -|- / d (/i) dx : |E| = 1 

i JPAElfi 

Hence, performing variations of 

Per(E) -k f d ih) dx 
JfaeI\ t-h 

within the class of sets of unit volume, for every 4* G C'^(M"'; M""), we again find 
(3.26), where X^^\t) is the Lagrange multiplier related to the constraint |E| = 
1. Observe that in both cases, we can choose a sequence of 'k G C'^(M”;M”) 
approximating e (h) on dE[^\ and conclude that X^^\t) is given by (3.27). □ 



It is now clear that the regular part of d*Ef^^ is of class Indeed, by 

choosing a suitable system of coordinates, d* e\^^ can be written in a neigh¬ 
bourhood of any regular point as the graph of a function u solving the 
following equation in the sense of distributions: 


div 


Vu 


= - X^^\ 


^Vl + |Vu|2^ 

Since is Lipschitz continuous, by standard elliptic regularity theory (cp., 
e.g., [18]), one then deduce that u G for every n G (0,1). 
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Next we prove that the whole family of sets defining the discrete flow up to 
time T > 0 is contained in a large ball, whose radius does not depend on the 
discrete time-step h but may depend on the T. 

3.4.3. Lemma (Boundedness of minimizing sets). Let be an ap¬ 
proximate solution and letT>0. Then there exist > 0 (depending on 

T, n, and Eq only) such that, if h < ho, then e[^'^ C for all t S [0, T]. 

Proof. We hx h > 0, and for every t G [0, T) we let 

rt := inf{r > 0 : e[^'* C Br}. 

We notice Br^ n dE[^'^ / 0 is made of regular points (because there are no 
singular minimizing cones contained in a half space, cp. [19, Theorem 15.5]), 
and moreover 

Br, n dEl'^^ C {y G dE\^'^ : F^(.)(y) > o} . 

By this observation and the Euler-Lagrange equation v^^\t,y) = X^^\t) — 
H (h)(y), it follows that 

rt<rt-h + h\\^'^\t)\. 


Iterating the above estimate, we then deduce that 

rT<'rQ+ [ \X^^\t)\dt VrG[0,r]. 


(3.28) 


To get some control on X^^\t) we consider T G Cg(M"';M”) such that 'I'(a:) = x 
in Br,, and, since for h small enough \E^^'^\ > ^ by (3.17), using T as test in 
(3.26), we obtain 


f lA'*!*)! < 


divT dx 


X^^\t) [ 

J E 

/ (diVr.Uh)'^ 

J^eW V dEl I 


X^^\t) 


IdE, 


(h) 


< (n — l)Per(E|^^) -|-rtPer(£'j^^)^/^||u^^^(t, • 






T dn 




n—1 


^dR 


n—1 




(3.29) 


where we used iTj < rt on dEj-^\ Integrating in time and using the Cauchy- 
Schwarz inequality together with (3.16) and (3.22), we obtain 

\ 1/2 


riA<'‘>(*)i 

Jo 


dt < Cn,0T -\- Cn,0 




(3.30) 


Combining (3.28) and (3.30), it follows that 

1/2 


Ft < I’D + CnfiT + Cn,o (^J rj dt^ for all r G [0, Tj. (3.31) 
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The remainder of the proof is a standard ODE argument. Indeed, squaring 
both sides of the equation and redefining Cn,o yields 

^ {tq -h r^) for all r G [0, T], 

where F{t) = f^r^dt. Integration in r over the interval [0,T] yields 

rfdt < Cn,0,T, 

and thus the statement in Lemma 3.4.3 follows via (3.31). □ 

3.4.4. Corollary. For every h > 0 small enough, it holds 

r f r\x^^\t)\^dt < cn,o,T. 

Jo JaEl^> Jo 

Proof. The integrability of follows from (3.29) and from Lemma 3.4.3; 
the one of H (h) follows taking into account the first variation (3.26) and the 
integrability of the discrete velocity, Lemma 3.4.1. □ 

For every /i > 0 we set 

S(h):={f: lEfl^l}. 

3.4.5. Corollary. For every h > 0 small enough, we have 

|S(h)| < Cn,0,Th. 

Proof. In view of Lemma 3.4.2, it is 

S(/i) C {f G [0,r) : |A('^)(f)l > 1/Vh} , 
and thus we have by Corollary 3.4.4 

|S(/i)| < h [ \X^^\t)\‘^dt < Cnfi,Th. □ 

Jo 

3.5. Proof of Theorem 2.3.1. The proof of the existence of a flat flow is 
now a simple consequence of the results above. Indeed, by (3.16), (3.17) and 
Lemma 3.4.3, one can find sets {Et}tgQ+, where Q+ denotes the set of positive 
rational numbers, and a subsequence | 0 such that 

lim AEt\ = 0 VfGQ+. 

fc^ + CX) 

Using the triangular inequality and Proposition 3.3.1, we deduce that 
\EtAEs\< lim 

fc^-l-oo V / 

< Cn,o\s - V0<s<fGQ+. (3.32) 

Now a simple continuity argument implies that the sequence converges to 
sets Et for all times f > 0 and satisfies (3.32) for all s, t G [0, -|-oo). Finally, note 
that passing to the limit in (3.15) yields that \Et\ = 1 and Per(i7i) < Per(i7o) 
(cf. [16, Sec. 5.2.1]) for a.e. t G (0,T). 
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3.6. Remark. It is also possible to show that the sequence of characteristic 
functions 


are precompact in L^((0,T) x M”) for every T > 0, thus giving an alternative 
proof of the theorem. 


4. Distributional solutions 

In this section we prove Theorem 2.3.2 on the existence of distributional 
solutions. The two main ingredients of the proof besides the estimates of the 
previous section are the hypothesis (2.2) on the continuity of the perimeters of 
the approximate solutions and the following proposition which links the discrete 
velocity to the distributional time derivative of the flat flow. 

4.0.1. Proposition. Let n < 7 and be an approximate solution to 

the volume-preserving mean-curvature flow. Then, for every f G +oo) x 

M”) it holds 


lim 

/n-o 


r+oo r r 

/ — d>dx — I (hdx 

Jh 2/l L J 


dt 



(fvhd'H^ ^dt 


= 0 . 

(4.1) 


Assuming the proposition we give a proof of the theorem. 


4.1. Proof of Theorem 2.3.2. It follows straightforwardly from (2.2) (cp., 
for instance, [3, Proposition 1.80]) that TC^~^\^d *weakly converges to 
for almost every t G [0,+oo). In particular this implies that, for 
almost every t G [0, +oo), the boundaries d *converge to d*Et in the sense 
of varifolds: namely, for a.e. t G [0, +oo) it holds 


lim 

k—^oo 


[ F(a 


.n^^h^){x)) dn 


n—1 


(x) 



VEt{x))dTr ^(x). 


(4.2) 


for every F G (70(1^” x M”’). Indeed, by a simple approximation argument it is 
easy to verify that it is enough to consider F G C'c(K"' x M"'). Then, for every 
e > 0 we pick a continuous function '■ such that 




dn^-^ 


< 
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and estimate as follows 


lim 

fc—)-+oo 


dn^ l(x) 


< lim ||Z)F||lo 
/ c^+co 


Jd-E\ 




n—1 


< lim ||DF||loc Per(4''"Y'^M f 

k^+oo" yjQ 

= WDFWlo. FeT{Et)y^ ( [ \UE, 
\Jd*Et 


1/2 


g(^fc) F-t 


EUhu) - k>F?‘ dW ^ 




1/2 


< \\DF\\l^ Pev{EtY^^s, 


where we used (2.2) and the weak convergence of the the vector valued measures 
V^ VEtdF'~^\-d*Ft in the following way: 

|2 j-i/n—1 


lim 


f 

'Jd*Ei 


k—^-\-oo Jf)* 


VUhk) - ^el dW 


/c—>- + CxD J Q* 

= [ + - 2vEt ■ Ve) d'H'^~^ 

Jd*Et 


2 Jn/n—1 


We/ - Ve\ dH 


Next we use Lemma 3.4.1 and Corollary 3.4.4 in conjunction with the results 
in Hutchinson [23, Theorem 4.4.2] to deduce the existence of functions v : 
[0, -|-oo) X M” —>■ M, A : [0, -|-oo) —M and H : [0, -|-oo) x M” —>■ M such that 

r |A|2 dt + r [ (jup + \H\^) n^-Ut < Cn,0,T, 

Jo J 0 J d* Et 


and 


/7 

Jo Jai 

fl 

Jo Jdl 


lim 

fc^ooJo JaE^^'^'* 

j-T 

lim 

Jo 


Vh^^fjidV.'^ ^dt = 


fl 

Jo Ja* 


vcjidT-F ^dt, 


Xih, 


'0 Jd*Et 

rT 


(4.3) 


^c/)dn^-^dt= f f Xc/)dn^-^dt, (4.4) 
Jo Ja*Et 

T 


lim f f H/H^)Ej,^H^)-^d'H^-^dt=f f H-^dn^-^dt, (4.5) 

k^°oJo JaE'^^k) F Jo Ja*Et 

for every cj) S (^^([OjT) x M”) and every $ £ (^^([O,T) x M"';M"'). 

In particular, testing (4.2) with F{x,v) := divT — v ■ VTi/ for some T G 
C^(]R'^;M”) and using (4.5), by a simple approximation argument we conclude 
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that, for a.e. t € [0,+oo), 

f div 

Jd*Et 


= lim 


1 


div. 


fc—>+00 




^dn 


n—1 


= lim 


1 


/c—>+00 J Q* 


(hk) 




ihk) 


dn 


n—1 


I 

Jd*Et 


VE, - ^HdU 


n—1 


thus showing that H{t, ■) is the generalized mean-curvature of Et for a.e. t € 
[0,-|-oo) (cp. (2.1)) and proving (2.3) of Theorem 2.3.2. 

Similarly, (2.4) and (2.6) follows from (3.26) and (3.27) by using (4.3) and 
(4.4). 

We need only to show (2.5). To this aim we use Proposition 4.0.1. For every 
(j) G ([0, -t- 00 ) X M”), by a change of variables we have that 


POO ^ r r 

/ / (hdx — Sdx 

Jh Ue^^^ Je^^\ ■ 


dt 



((/>(t, x) — (f>(t + h, x)) dx dt 


-h 



where we used that = Eq for t G [0,/i). Therefore it follows by a simple 
convergence argument that 


lim 

h^O 


I \ \ f (pdx- [ 

Jh h 1 Je^^'> Je 


(h) 


(pdx 


dt 



-^(t, x) dx dt 
dt^ ' ’ 



In view of (4.3) and (4.1), we conclude (2.5) straightforwardly. 


4.2. Tilting of the tangent planes. In this subsection and in the next one 
we give the proof of Proposition 4.0.1. We follow closely the arguments in [26] 
and for the sake of completeness we provide a detailed proof in different steps. 

This subsection is devoted to the estimate of the tilting of the normals around 
points of small curvature. We recall that we assume in this section n <7 (in 
particular, the approximate solutions of the volume-preserving mean-curvature 
flow are everywhere of class C'^’'^). 

4.2.1. Lemma. For given constants ^ < (5 < a < 1, there exists a con¬ 
tinuous increasing function cj : [0,1] —)• M with ti;(0) = 0 with the following 
property. Let t G [2h, -|-oo), {El.^^}t>o be an approximate solution to the volume¬ 
preserving mean-curvature flow, and let xq G dE[^^ he such that 

^2/ei?^^^(a;o)n(i?f AE^). 


v^^Ht,y)\<h--^ 


(4.6) 
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Then there exists v £ M” such that |z/| = 1 and 

WgEwiy) -n\< uj{h) V ?/ £ Bj,p{xq) n dEf‘\ (4.7) 

(y) -^\< yyGBf^p (xo) n dEi^{. (4.8) 


Proof. Let 0 < R < h 2 ^ and let F C ME be any set such that E^'^IS.F CC 
BRhi^{xo)- By the minimizing property of we have that 

Per(£’f ^, BehP (a^o)) < Per(F, B^hP (a^o)) + v / , d (i.) {y) dy 

h JFAEi^^ t-h 


+ ^(||F|-i|-||p;f)|-i|). 


(4.9) 


A straightforward computation yields 


||F| - l| - I IE) '’I - l| < |EAE^ 


Wi 


1 

h 




dpW iy)dy < 


7n + 1 
y/h 


FAE, 


Wi 


where we used that \v^^\t,y)\ < Rh^ ^ + 7 n h V 2 < (-y^ + 1 ) h 2 for all y £ 
BRhP{xQ)f^{Ef^'^ AF) thanks to the fact that xq £ dEf^\ Proposition 3.2.1 and 
the 1-Lipschitz continuity of the signed distance sd„(i,) . Combining the above 

^t-h 

estimates with (4.9), we obtain 

Per(Ef ,E^,,(xo)) < Per(E,E^,,(xo)) + |FAEf)|. 


y/h 


(4.10) 


Next we introduce the sets 



[z £ : 

y-xo 
^ hP ’ 

y £ Ef)} 

X.IY ■■= \ 

[z £ : 

y-xo 

^ hh ’ 

y e e£1} 


By a simple rescaling argument, from (4.10) and from the analogous estimates 
at time t — h (recall that t > 2h) we deduce that for s = t,t — h 

Per(Ef Er) < Per(E, Br) + ( 7 ^ + 2) |FAEf 

VEAE^’^CCEr. 

This implies that and are both (A/^, r/j)-minimizers of the perime¬ 

ter on A/j := ( 7 ^ -|- 2 )h ^~2 and := /i 2 “^. By the precompactness for se¬ 
quences of A/j-minimizers (cf. [27, Prop. 21.13]), we conclude that we can find 
a subsequence (not relabeled) verifying 

lim XeW,p = XeP and lim x (/^),/3 = Xr/? in ^ioc(I^’^)- 

Moreover, using the lower semicontinuity of the perimeter with respect to 
convergence and /3 > ^, we deduce that E^ and E 2 are locally minimizing 
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the perimeter. By the assumption n < 7 and a Bernstein theorem (see [19, 
Theorem 17.3]), E^, are half-spaces. Moreover, by hypothesis it holds 

thus implying that E^ = E^, and by the fact that both are hyperplanes there 
exists € M” with |i/| = 1 such that 

E^ = E^ = {z£R^ : z-u<0}. 

To reach the conclusion of the lemma we need only to invoke the regularity 
theory of A-minimizing set (cp. [27, Theorem 26.3]) and conclude that dEg^^’^ 
is uniformly C^’'^ in Bi for s = t,t — h, thus leading straightforwardly to (4.7) 
and (4.8). □ 

4.2.2. Corollary. Under the hypotheses of Lemma let C^/3/2(xo, z^) C 

M”' be the open cylinder defined as 

:= jx G M"" : [(a; — xq) ' i/|x — xqP — |(x — xq) ■ z/|2 < ^'’/ 2 | . 

Then, there exists a dimensional constant C > 0 such that 


/ {Xp(h) - XpW )dx- 

ic,,/2(xo,p) Ja 


sd (h) dTL 


n—1 


< Cuj 


(h) [ 

Jc 




\Xp(h)-Xpw\dx. (4.11) 


Proof. From Lemma 4.2.1 we know that, for h sufficiently small, dE[^'^ 
and dE\% in i 2 {xq,v) can both be written as graphs of functions of class 
Namely, by an affine change of coordinates we can assume without loss of 
generality that xq = 0 and v = Cn, and for simplicity we set C := C ^ 0 / 2(0 Tn)- 
With this assumption we then have that for s = t,t — h 


n C = {{yjsiy)) G X M : |y| < hy 2 } , 
where fs : C —)• M are C^’'^ functions with 

fs\\L°°{B,^l3/2) - 

In view of Fubini’s theorem it is then clear that 



\My) 


ft-h{y)) dy, 


ft-h{y)\ dy. 


(4.12) 
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Moreover, from (4.12) it follows that there exists a geometric constant C > 0 
such that, for every y G 

|sd^(M (y,/t(y))Vl + \Vft{y)? - ift{y) - ft-h{y))\ < Cu:{h) \ft{y) - ft-h{y)V 
Therefore, one infers (4.11) as follows 


L 


(h) 

dE'>'^>r\C 


sd^rw d'H” ^ — 


IB 


iftiy) - ft-h{y))dy 


hf>h 


[ (h) (y, ft{y))^/TTWMyW - (ftiy) - ft-h 


dy 


< Coj 

where we used that (4.12) 


(h) f \ft-ft-h\dy, 


□ 


We are finally ready for the proof of Proposition 4.0.1. 


4.3. Proof of Proposition 4.0.1. We fix any time t G [2h,+oo). For every 
xo G dE[^\ we fix a G (^, 2 {r^i) ) consider the following open set A^q 
defined as follows: 

(i) if (4.6) holds, then we set := C^/9/2(xo, i^) where u G M” is the unit 
vector in Lemma 4.2.1; 

(ii) otherwiese we set A^^ := 

Note that by Proposition 3.2.1 we have that is a covering 

of Moreover, by a simple application of Besicovitch’s covering 

theorem, cp. [16, Ch. 1.5.2] (applied, for example, to the balls to B^/3/2(xo) C 

Axq), there exists a finite collections of points I C dE[^^ such that {^xqIxoS/ 
is a covering of E^^'^ AE^^^. 

We estimate the contribution of the integrals in (4.1) in every A^q with xq G / 
in two steps, depending on whether (i) above applies or (ii). 




sd (h) d'H" ^ 

^t-h 


Estimate in case (i). We use Corollary 4.2.2 and deduce that 
/ {Xp,(h)-Xp,(h))4>dx- sd(h)4>d'H'^~^ 

< l</'(a;o)| / {Xp(h)-Xp(h))dx- 

Ja,^ Jd 

+ / {XpW - Xpw) i4>-4>{xo))dx- 

JA^^ Jd 

< C {uj{h) ||(/>||loo -F hd ||V())||loo) f \x (h) - V (h) I (4.13) 

/ A -^t — h 

where we used the fact that Axq = C^Pj 2 {xo, v). 




sd (h) i4> - 4>{xo)) dH 
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Estimate in case (ii). By assumption there exists a point ^i^o) ^ 

such that \v^’^\t,yo)\ > h°‘~^. Without loss of generality we can 
assume that yo € E^^'^ (the other case can be treated analogously and we leave 
the details to the reader). It is then clear that Bh‘^/ 2 {yo) C M” \ E^^f^ and 
v^^\t,y) > for every y £ Bh°‘/ 2 {yo)- Since /i“ < 27 ^/ 1 ^/^, we can apply the 

density estimate in (3.11) and deduce that 

< / \v’^\dx. (4.14) 


Similarly, by the density estimate in (3.12) and Proposition 3.2.1 we deduce 
that 


IB 




(h) 


d^^ <Cnh"/\ 


J. ■T'h 
^t-h 


(4.15) 


From (4.14), (4.15) and Bh°‘/ 2 (,yo) C we then deduce that 



IXg('i) Xr('i) 


+ 


d (h) dH 


n—1 






L 




Iv^ldx. 


(4.16) 


We can then sum (4.13) and (4.16) over xq G I and, recalling (3.15), (3.16) 
and (3.17), we get 


1 (x^(« - x^(.)j 4‘d^ - Bd^»Jdn 

sE|/^ j. 


sd (/i) (fidJi 


n—1 




xoei •'^^0 

< CnflUh) Uh^+hf^ ||V(()|Uco + r/^-(-+i)“+i Wd^u. 


|£;f )Ai?W I + 


L 


\v^\ dx 


Ie[’^^ae[^X 

where we used the finite finiteness of the covering. 


(4.17) 


Finally, integrating in time and using (3.18) and (3.19) we get 

r+00 1 r r r 1 ^+00 p 

/ h dt- / 4>Vhdn'^~^dt 

J2h h 1 E JdE[^ 

< C'„,o,rfw(/i) UWl^ + IIV^IUoc + h"/^-in+l)a+l II^II^ 


where T > 0 is such that supp((/>) C [0,r] x M"'. Recalling the definition of ol 
and taking the limit as h goes to 0, we conclude (4.1). 










Volume-preserving Mean-curvature Flow 


23 


Appendix A. A density lemma 

We premise the following density estimate for one-sided minimizers of the 
perimeter. The estimate can be easily deduced from the original arguments by 
De Giorgi exploited for minimizers [14]. 

A.0.1. Lemma. There exists a dimensional eonstant Cn > 0 with this prop¬ 
erty. Let E C M”' be a set of finite perimeter, R, n > 0 and xq € dE he such 
that 

Fer{E) <Fer{E D Br{xo))n\Br{xo) \ E\ V 0 < r < i?. (A.l) 

Then, 

Cnr'^ <\Br{xQ)\E\ V 0 < r < min {i?,. (A.2) 

Recall that xq G dE if min{|ilr(3:o) \ £^|, n ilr.(3:o)|} > 0 for every r > 0. 

Proof. Without loss of generality, we can assume xq = 0. We use the 
following indentity which are true a.e. r > 0: 

FeviEUBr) =n^-^{dBr\E) + Fer{E,R^\Br{x)), (A.3) 

Per(R^ \E)= n^-^{dBr \ E) + Per(E, Br), (A.4) 

Per(p;) = Fer{E,Br) + Per(£;,M’" \ Br). (A.5) 

Indeed, if E were smooth, these formulas follow for all the r such that Br and 
E have transversal intersections. Otherwise one can argue by approximation. 
Using now (A.l), we deduce that, for a.e. r > 0, 

Per(R^ \ E) ^ R^-^dBr \ E) + Fei{E, Br) 

(A.5) 

< n^-^idBr \E) + Fer{E) - Fer{E, \ Br) 

(A.l) & (A.3) 

< 2n^-\dBr\E)+ n\Br\E\. (A.6) 

By the isoperimetric inequality [19, Corollary 1.29], there exists a dimensional 
constant G > 0, such that 

n-l (A.6) 

C\Br\E\— <Fer{Br\E) < 2'H^-^{dBr\ E) + n\Br\ E\. (A.7) 

Setting /(r) := \Br \ ^l, by the coarea formula [16, 3.4.4], it holds 
TE~^{dBr \E) = f'{r) for a.e. r > 0. 

Hence, (A.7) reads as 

C f{r)~ <2f{r) + nf{r). (A.8) 

Finally, note that /(r) < ojnr'^, from which /(r) < r f{rY~^^'^. Therefore, 
there exists a dimensional constant Cn > 0 such that if 0 < r < min{i7, Cn 
then the last term in (A.8) can be absorbed in the left hand side and deduce 
that 

Tl 1 

/(r)— <C/'(r). 
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Integrating (A.7) (recall that /(r) > 0 for every r > 0) we get the desired (A.2) 
for every 0 < r < min{i?, and, by changing the dimensional constant 

Cn > 0, for every 0 < r < min{i?,□ 
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